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j^^ 1. Introduction 

^7^ The nuclear symmetry energy, Eg — 2^52", where S — (n„ — np)/{n = rin + rip) 

^ with n„ and Up denoting the neutron and proton densities and E = E{n, 5) is the 

[v>< energy per particle, measures the stiffness encountered in mafcing a system of nu- 

^ cleons isospin-asymmetric. Figure [T] schematically shows how the symmetry energy 

connects several nuclear and astrophysical observables. The difference between the 
energy per baryon of pure neutron matter and that of symmetric nuclear matter 
(containing equal numbers of neutrons and protons) at any particular density is 
largely given by the density dependent symmetry energy. Below, we highlight some 
recent work that has shed light on some of the connections between the symmetry 
energy and data from terrestrial experiments and astrophysical observations. 

2. The Skin Thicknesses of Heavy Nuclei 

Traditionally, constraints on the nuclear symmetry energy have been derived from 
mass measurements of nuclei. For neutron star physics, the correlation that exists^''^ 
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Fig. 1. The nuclear physics observables (top panels) and astrophysical observables (lower panels) 
are both connected to the nuclear symmetry energy, which is determined from nuclear many-body 
physics (adapted from Ref. 1). 

between the neutron skin thicknesses in heavy nuclei, SR (the difference between 
the neutron and proton root- mean-square radii) , and the pressure P of pure neutron 
matter at a density of n ~ 0.1 fm~'^ is particularly useful. Accurate measurements 
of 5R can establish an empirical calibration point for the pressure of neutron star 
matter at subnuclear densities. The connection between the neutron skin thickness 
and the symmetry energy has been known from Bodmer's work in the 60's.* The 
Typel-Brown correlation between SR and P(5no/8), which is closely related to the 
density derivative of the symmetry energy, demonstrates clearly the new information 
that could be obtained by accurate measurements of skin thicknesses in heavy nuclei. 
Figure [2] shows the correlation between the skin thickness SR of ^°*Pb and the 
pressure of beta-equilibrated matter at 0.1 fm^"^ for several potential models (based 
on the Skyrme interaction) and field-theoretical models. "'^ The skin thickness of 
208p|^ is scheduled to be measured at the Jefferson Lab in the summer of 2008 in 
the PREX experiment'^'^ and will likely provide a stringent constraint. 



3. Intermediate-Energy Heavy-Ion Collisions 

The possibility of determining the equation of state (EOS) of nucleonic matter from 
heavy- ion collisions has been discussed for almost 30 years. A number of heavy-ion 
collision probes of the symmetry energy have been proposed including isospin frac- 
tionation,^'^ isoscaling,^' ^'^ neutron-proton differential collective flow,^-'^ pion pro- 
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Fig. 2. The correlation between neutron skin tliickness of 208pjj ^^^^ ^j^g pressure of neutron star 
matter for the models described in Ref. 1. 



duction/^ isospin diffusion/'^ and neutron-proton correlation functions.^* Determi- 
nation of the EOS from heavy-ion data involves comparisons of transport model 
simulations^^^^^ with experimental data given an input EOS. Transport model sim- 
ulations track the evolution of the phase space distribution function (not necessarily 
that of an equilibrium distribution function at finite temperature) at momenta that 
usually exceed those found in the static initial configurations. As the driving force is 
the density functional derivative of the energy density (the in- medium cross sections 
control the collision integral), access to the cold EOS at high densities is afforded. 
Isospin diffusion is caused by the exchange of neutrons and protons between nu- 
clei in a heavy-ion collision (Fig. pi). This diffusion process, driven by the symmetry 
energy, moves the target and projectile nuclei toward isospin symmetry. To gain 
access to details of the diffusion process, fragment emission during and after the 
collision must be taken into account. This requirement is achieved by considering 
the ratio^* 

_ 25^+B - S^+^ - S^+B 

where A and B denote nuclei with different isospin asymmetries and S is the isospin 
asymmetry of the projectile- like fragment. 

Recently, isospin diffusion has been exploited to constrain the symmetry en- 
ergy from reactions involving ^^^Sn and ^^''Sn at the NSCL,^-^'^*^ leading to Rs ^ 
0.46. In conjunction with an isospin- and momentum-dependent transport model, 
IBUU04,^^'^^ the NSCL data on isospin diffusion can be used to constrain the sym- 
metry energy. Writing the symmetry energy in terms of a kinetic part and a potential 
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Fig. 3. Isospin diffusion process during a heavy-ion collision (from Rcf. 19). 



part as 



Esin) = S„{nlnoV , (2) 

where S.^ is the symmetry energy at the nuclear equihbrium density, tiq = 0.16 fm^ , 
the constraint 0.69 < 7 < 1.05 has been found. ^^'^'^ This constraint is consistent with 
the symmetry energy inherent in the EOS (computed using Monte Carlo simulations 
with input two- and three-body interactions which are matched to nucleon-nucleon 
scattering phase shifts and the energy levels of light nuclei) of Akmal, et al.-^"* 
(APR). This constraint also rules out models with values of 7 > 1.05 found in some 
field-theoretical models. 

Because intermediate-energy heavy-ion collisions provide a constraint on the 
symmetry energy at the same densities as would be probed by a measurement of 
the neutron skin thickness of ^°^Pb, the NSCL data also provides a restrictive range 
for its neutron skin thickness. This connection was used in Refs. 25 and 23 to show 
that the neutron skin thickness of ^"*Pb should be at least greater than 0.15 fm in 
order to be consistent with the NSCL data, with values between 0.23 fm and 0.27 
fm favored by the transport model simulations (Fig.Q. 



4. Neutron Star Radii 

Neutron star radii tend to probe the density dependence of the symmetry energy 
around the nuclear equilibrium density, rig = 0.16 fm^ . Lattimer and Prakash^^ 
found that the radius i? of a neutron star exhibits the power law correlation: 



i?~C(n,A/) [P{n)f 



(3) 



where P{n) is the total pressure inclusive of leptonic contributions evaluated at 
a density n in the range tiq to 2no, and C{n,M) is a number that depends on 
the density n at which the pressure is evaluated and on the stellar mass M. The 
left panel in Fig. 5 shows this correlation as RP^" versus R for stars of mass 
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Fig. 4. The correlation between the neutron skin thickness and the amount of isospin diffusion 
in coUisions between Sn isotopes as measured at the NSCL (from Ref. 23). 



1.4 M0. Neutron star radius measurements, especially those with uncertainties less 
than about 0.5 km, constrain the symmetry energy above the nuclear equilibrium 
density. These constraints will be much improved when simultaneous mass and 
radius measurements of the same object become available. 

If there is no phase transition between n^ and a few times no, the range in 
which neutron star radii are determined mainly by the symmetry energy, results 
from the isospin diffusion data at the NSCL can be used to constrain neutron 
star radii. ^^ As only EOSs with symmetry energies between a; = and x = —1 
(where a; is a parameter designed to vary the density dependence of the symmetry 
energy without modifying the magnitude of the symmetry energy at no or the 
isospin-symmetric part of the EOS) are consistent with the isospin diffusion data, 
this range of x values is representative of the possible variation in neutron star 
structure that is consistent with terrestrial data. Neutron star radii, while being 
strong functions of the symmetry energy, are also affected by contributions from the 
isospin-symmetric part of the EOS, especially at high densities. About 5% difference 
is representative of the radius uncertainty stemming from the symmetric part of 
the EOS. The conclusion of Ref. 27 is that only radii between 11.5 and 13.6 km (or 
radiation radii between 14.4 and 16.3 km) are consistent with the a; = and x — —\ 
EOSs (see the right panel in Fig. 5). 
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Fig. 5. (Left) The value RP~°' for several field-theoretical and potential models considered in 
Ref. 1. The values of a correspond to the density indicated in the upper right corner of each block. 
(Right) The constraint on the radius of 1.4 solar mass neutron stars as determined from isospin 
diffusion in heavy-ion collisions (from Ref. 27). 



5. The Direct Urea Process 

The long-term cooling of a neutron star is chiefly determined by its composition. 
Beta equilibrium and charge neutrality determine the proton fraction in neutron- 
star matter and thus the critical density for the onset of the direct Urea processes 
n —^ p + e + D and e + p —^ n + v, which cool the star more rapidly than the 
modified Urea processes in which an additional nucleon is present. ^^ The direct 
Urea processes, however, require a sufficient amount of protons in matter (of order 
10-14%). 

Larger symmetry energies induce larger proton fractions (with matter being 
closer to isospin-symmetric) and smaller critical densities for the onset of the direct 
Urea processes. However, higher than quadratic terms in the energy E{n, S) of 
isospin asymmetric matter can have an important role to play.^^ Recently, Steiner^° 
has shown that quartic terms in E{n, 6) play an important role in determining the 
critical density for the direct Urea process. Such terms can be easily generated 
within the context of field-theoretical models. This is demonstrated in Fig.[6]for the 
Akmal, et al., (1998) EOS, for which the relative size of the quartic and quadratic 
terms is parametrized by i] as described in Ref. 30. The mass and radius are virtually 
unchanged, whereas the threshold density for the direct Urea process changes by 
more than a factor of two. 

Gusakov, et al.,^-*^ have investigated the cooling of neutron stars using the EOS 
of APR.^"' They found that, because of the direct Urea process, stars with masses 



February 5, 2008 4:33 



WSPC - Proceedings Trim Size: 9.75in x 6.5in proc3 



larger than about 1.7 M0 cool so rapidly as to be cooler than nearly all of the 
observed neutron stars. Our work offers a possible resolution: quartic terms can 
play a role at high density to turn off the direct Urea process thus making the 
computed cooling curves match the comparatively warm neutron stars. 
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Fig. 6. The maximum mass, tlie radius of the maximum mass star, and the critical density for 
the direct Urea process as a function of rj, which describes the strength of quartic terms in the 
symmetry energy (from Ref. 30). 



6. Outlook 

In addition to the role of the symmetry energy in the few areas highlighted here, 
its importance in controlling the cooling times of transient x-ray bursters, seismic 
activity of neutron star surfaces, ejection of baryons during binary mergers, etc., is 
only beginning to be appreciated. The PREX experiment, heavy-ion experiments, 
neutron star mass, radius and surface temperature measurements, observations of 
transient x-ray bursters etc., all hold keys to pin down the magnitude and density 
dependence of the symmetry energy in addition to delineating its pervasive role. 
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